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Abstract 
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1 Introduction 

The inclusive production of high-transverse-momen- 
tum single hadrons by collisions of protons and the 
quasi-real photons emitted by an electron beam is an 
interesting process for testing the QCD-improved par- 
ton model. First experimental results on this process 
have been presented recently by the HI Collabora- 
tion at HERA [1]. These data have been compared 
already with theoretical predictions to leading order 
(LO) in QCD. The LO QCD formalism consists of us- 
ing the tree-level results for the hard-scattering cross 
sections, the one-loop expression for the running cou- 
pling constant, and structure functions and fragmen- 
tation functions generated by one-loop evolution ker- 
nels. It is well known that such LO estimates suffer 
from normalization uncertainties due to the renormal- 
ization/factorization scale dependence. The normal- 
ization ambiguities are substantially reduced in the 
next-to-leading order (NLO) formalism. 

There exist two distinct mechanisms which con- 
tribute to the inclusive photoproduction of hadrons at 
high energies. The photon can interact either directly 
with the partons originating from the proton, giving 
rise to the "direct" process, or via its quark and gluon 
content, the so-called "resolved" process. In both cases 
high-px final partons are produced, which fragment 
into single hadrons with the same pt. 

Both contributions are of the same order in the 
perturbative expansion. The LO contribution of the 
resolved part is given by LO parton-parton scatter- 
ing terms, multiplied by LO structure functions of the 
photon and proton and LO fragmentation function of 
the final hadron. Since the parton-parton scattering 
terms are of O (a 2 ) and the photon structure function 
of 0(a/a s ), both LO contributions of the direct and 
resolved parts are of 0(aa s ). 

Recent data from the HI and ZEUS Collaborations 
[1,2] indicate that the resolved photoproduction domi- 
nates at HERA energies, for small transverse momenta, 
Pt , and negative rapidities of the produced hadron rel- 
ative to the laboratory system, yi a b (we take yi a b to 
be positive in the direction of the electron beam). So 
either by appropriate kinematical cuts in p^ and yi a b 
or by observing the remnants of the resolved photon 
near the beam pipe it should be feasible to obtain in- 
formation on the resolved part alone [1]. Here it is 
understood that pt is still outside the region of soft 
particle production, which dominates for pr 1 GcV. 
Thus considering resolved photoproduction of single 
hadrons as a separate process, it is important to obtain 
quantitative predictions for it. For this purpose, we 



have calculated NLO corrections, which are of O(a^) 
for the parton-parton scattering process, in the case of 
the resolved photoproduction alone. It is clear that in 
the high-px range direct and resolved photoproduction 
yield comparable contributions, which must be com- 
bined in order to obtain a reliable prediction. This will 
be left for future work. 

To perform a consistent calculation of inclusive 
single-particle production in the NLO formalism, we 
need the hard-scattering cross section in NLO, with 
two- loop a s (/it), two- loop evolved structure functions of 
the photon and the proton, and two-loop-evolved frag- 
mentation functions. Only in conjunction with each 
other are these elements unambiguously defined. 

As far as the proton structure functions are con- 
cerned, NLO parametrizations exist already for some 
time [3,4,5]. As for the photon structure functions, 
recently three NLO sets of parametrizations have ap- 
peared in the literature [6,7,8]. They have been tested 
against the available experimental data on the photon 
structure function ^(x, Q 2 ). 

Unfortunately, such analyses do not exist for the 
fragmentation functions of quarks and gluons in gen- 
eral. Only the 7r° fragmentation functions in NLO have 
been developed recently [9]. This has not been done 
for charged particles yet, in which we are mainly in- 
terested here. The existing information on fragmenta- 
tion functions of charged particles is based on fits to 
data from e+e~ annihilation and deep-inelastic lepton- 
nucleon scattering in connection with LO evolution 
[10,11]. These LO fragmentation functions have been 
tested in inclusive charged-particle production in pp 
scattering in a previous communication [12], henceforth 
referred to as I. Reasonably good agreement with ex- 
isting data from the UA2 and CDF Collaborations has 
been found. Since the center-of-mass (CM) energies 
in these experiments, yfs = (540-1800) GeV, are even 
higher than the HERA CM energy, y/s = 298 GeV, we 
feel confident that these LO functions give a reasonable 
account of the fragmentation of quarks and gluons into 
charged hadrons. 

The NLO corrections to the parton-level cross sec- 
tions appropriate for inclusive single-particle produc- 
tion have been worked out by Aversa et al. [13]. They 
provide the basis of the NLO corrections to charged- 
particle and neutral-pion production presented here. 

In this paper we shall predict cross sections of 
single-charged hadron and single- ir° production in low- 
Q 2 ep reactions at the HERA energy. This work is an 
extension of I, where inclusive hadron production in pp 
collisions was studied. The formalism used here is quite 
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similar. One has to adjust the initial state and take 
care of the fact that the incoming electron produces 
a spectrum of photons, which is usually described in 
the Weizsacker- Williams approximation (WWA) [14]. 
Compared to pp collisions, we expect in jp processes a 
different behaviour of the cross section due to the pres- 
ence of the point-like component in the photon struc- 
ture function. 

Experimental data of charged-particlc and ir° pro- 
duction at relatively low CM energies come from fixed- 
target experiments at CERN [15]. In these experiments 
the resolved-photon component could not be separated 
from the direct one. Resolved photoproduction cross 
sections for single hadrons have already been measured 
at HERA [1]. A better statistic is now being accumu- 
lated and a wider range of pt is under study. These 
data will probe the QCD-improved parton model in an 
energy regime that lies between y/s = 200 GeV, the pp 
CM energy of UAl, and y/s = 540 or 630 GeV as used 
by UAl, UA2, and CDF. 

The outline of this paper is as follows. In Sect. 2, 
we give a brief introduction to the theoretical input 
and the structure functions used for the calculation. 
Numerical results are presented in Sect. 3, where we 
compare with the first data from the HI Collaboration. 
Section 4 is reserved for a discussion of the results and 
some outlook to future work. 



2 Theoretical ingredients 
structure functions 



and 



actual electron acceptance must be incorporated. This 
will affect our results only through a normalization fac- 
tor of (2.1). Then in the WWA the inclusive single- 
particle cross section is related to the corresponding 
photoproduction cross section in the following way: 



E, 



d 3 cr(ep -» hX) 
d 3 p h 



f 



dxf l/e {x) E h 



d 3 cr(-/p -> hX) 
d 3 p h 



(2.2) 

On the right-hand side of (2.2), xE e is substituted for 
the photon momentum and x m ; n is expressed in terms 
of transverse momentum, px, and CM rapidity, j/cm, 
of the produced hadron as 
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The rapidity measured in the laboratory frame, yi a t>) 
is related to ycM by 



yiab = ycu - .in — 



(2.4) 



where the energies of the electron and the proton, E e 
and E p , are in the laboratory frame. According to 
current HERA conditions, E e — 27 GeV and E p = 
820 GeV, so that y\ a b = ycu — 1.71. 

In the QCD-improved parton model the single- 
particle photoproduction cross section on the right- 
hand side of (2.2) is expressed in LO as a convolu- 
tion of the LO parton-parton scattering cross section 
with the scale-dependent structure and fragmentation 
functions. We fix the momenta of photon, proton and 
outgoing hadron (h) by 



When we speak of photoproduction with almost real 
photons we have in mind the class of ep processes where 
the electron is a source of nearly massless, collinear 
photons, which collide with the opposing proton beam. 
The energy spectrum for the resulting photons is cal- 
culated in the WWA. In this approximation the longi- 
tudinal- and transverse-momentum components of the 
photon distribution decouple. The transverse momen- 
tum is integrated out with an angular cut of 6* max = 5° . 
Then the distribution in the longitudinal-momentum 
fraction, x, of the outgoing photon has in the leading- 
logarithm approximation the following form [14] : 



, , s a 1 + (1 - x) 2 E e 6 max 

f l/e {x) = In 

7r x m e 



(2.1) 



where E e is the electron energy, m e is the electron 
mass, and xE e is the photon momentum. Equation 
(2.1) determines the photon luminosity used in our cal- 
culation. When comparing with experimental data, the 
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(2.5) 



In this way, the notation is similar to the one used in 
I. Mutatis mutandis, the LO cross section is given by 
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where fcj = x 1 p 1 , kj = x p p p , and kh — Ph/xh are the 
parton momenta. The indices i, j, I run over gluons and 
Nf flavours of quarks. We assume Nf = 4 throughout 
our calculation and include the charm-quark thresh- 
old as implemented in the photon structure function 
parametrizations used. Fj(x p , M 2 ) is the proton struc- 
ture function appropriate for parton j and D^(xh, M 2 ) 
is the fragmentation function of parton I into hadron 
h. M is the factorization scale of the photon and pro- 
ton structure functions, which we assume to be equal, 
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Mf is the scale of the fragmentation function. Finally, 
d 3 a® k ,_ tk describes the process i + j — > I + X in 

0(a 2 (>))- 

The calculation of the NLO inclusive cross section 
is described in I. The relevant formulae are used with 
the replacement F? 1 (x\, M 2 ) -> F^(x 1 ,M 2 ) and 
F^{x 2 ,M 2 ) Ff(x p ,M 2 ). For the full calculation 
we need the structure functions of the photon and the 
proton and the fragmentation functions in NLO. As 
already mentioned in the introduction, fragmentation 
functions exist only in LO, except for the recently con- 
structed 7T° set [9]. We shall use the same fragmenta- 
tion functions as in I, namely those of Baier et al. [10] 
and Anselmino et al. [11] with the constraints spec- 
ified in I. Since these fragmentation functions corre- 
spond approximately to the MS factorization, we have 
also selected the MS scheme for the proton and photon 
structure functions. 

The proton structure function parametrizations are 
again taken from the package by Charchula [17]. For 
most part of our analysis we use the Morfin-Tung (MT) 
set Bl in the MS version. 

For the NLO photon structure functions we have 
three sets at our disposal [6,7,8]. All three have been 
tested against the available data on the photon struc- 
ture function T^{x, Q 2 ). Unfortunately, this leaves still 
appreciable freedom, in particular with respect to the 
gluon content of the photon. For our main results we 
use the sets by Gliick, Reya, and Vogt (GRV) [6] with 
Np = 4. In their analysis the Q 2 evolution starts with 
a valence- like input at a rather small scale, Q 2 ,, which is 
Ql = 0.25 GcV 2 in the LO fit and Q 2 = 0.3 GcV 2 in the 
NLO fit, respectively. Their NLO set refers to the so- 
called deep-inelastic 7 scattering (DIS7) scheme. For 
consistency in our analysis, we convert the GRV DIS7 
set into the more familiar MS scheme. For compara- 
tive studies at the LO level of the hard scattering we 
use also the LO set by GRV. Besides, we make use 
also of the LO set by Drees and Grassie (DG) [18] and 
the NLO set by Aurenchc et al. (ACFGP) [8] to inves- 
tigate the sensitivity of the cross section to different 
choices. The boundary conditions used by ACFGP are 
similar to those of GRV and the evolution is started 
at Ql = 0.25 GeV 2 . The gluon content of the pho- 
ton is particularly low in the ACFGP set. The charm 
distribution is totally described by the lowest-order ex- 
pression, which has its starting point at Q 2 , = 2 GeV 2 . 
Gordon and Storrow (GS) [7], on the other hand, start 
the Q 2 evolution at Q 2 , = 5.3 GeV 2 , which is consider- 
ably larger than the values of GRV and ACFGP. 

In the next section we shall present the results in 
two steps. In the first step, we assume that D^(xh, M 2 ) 



= 5(1 — Xh) for all hadrons h and partons I. In this 
way, we study the cross sections for the production of 
single partons. In LO they are identical to the single- 
jet production cross sections. In NLO they depend on 
the factorization scale M 2 and cannot be compared to 
jet production cross sections, since the recombination 
algorithm for defining jets is not built in. However, the 
cross section with <5-function fragmentation are simpler 
to calculate and reveal several features that are also 
relevant for the complete single-particle cross sections 
considered in the second step. 



3 Numerical results 

3.1 Single-parton inclusive cross 
sections 

In this sub-section we approximate all fragmentation 
functions by <5 functions. This leads to single-parton 
inclusive cross sections. For a first orientation, we dis- 
entangle the various parton-parton channels that con- 
tribute in LO to the hard-scattering cross sections or- 
dering them according to their initial state, viz. 
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where q and q' denote different quark flavours. The 
first parton originates from the photon, the second 
from the proton. These four classes of channels are also 
compared with direct photoproduction in LO, which 
proceeds through 75 — > qq and 7(7 — * gq, qg. The 
comparison is done for d 2 a/dydpT at y/s = 298 GeV 
(this corresponds to 27 GeV electrons colliding with 
820 GeV protons according to present HERA condi- 
tions) and yiab = —2 (with yi a b > in the direction of 
the electron momentum). For this calculation we use 
the LO set SL of the MT proton structure functions 
and the LO set of the GRV photon structure functions. 
We also fix the two scales \i and M to be \i = M = px- 

In Fig. 1 we show the px dependence of the chan- 
nels l)-4) contributing to the resolved -fp cross section 
and the px dependence of the direct process. We ob- 
serve that at small pt channels 1), 2), and 3) dominate 
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20 30 
Pt [GeV] 

Fig. 1. Transverse-momentum distributions of the various 
channels for resolved photoproduction and of direct photo- 
production in LO at y/s = 298 GeV and yi a b = —2 evalu- 
ated with the LO sets by GRV and MT-SL. The resolved- 
photon channels are grouped according to their parton- 
parton initial state; the first parton originates from the 
photon and the second from the proton 

by one order of magnitude over class 4) and the di- 
rect mechanism, whereas at larger pt (pt ~ 28 GeV) 
class 4), in particular the channels involving quarks 
only, contribute significantly. It is exactly in this pt 
range, more precisely for pt > 35 GeV, that the direct 
process exceeds the sum of the resolved processes. The 
pattern observed for the various resolved photoproduc- 
tion channels is similar to the situation encountered in 
I in the case of pp collisions: the gluon components 
of the photon and proton structure functions influence 
strongly the cross section at smaller pt, whereas the 
pure quark channels dominate at larger pt • 

In Fig. 2, the rapidity distributions of the above 
channels are displayed for fixed p T — 5 GeV. For 
yiab < 0, i.e., in the proton direction, classes 1) and 3) 
dominate and the direct mechanism comes in only in 
the extreme forward direction, at yi a b ~ 1. Of course, 
the relative size of the various classes in Figs. 1,2 
may change when realistic fragmentation functions are 
taken into account. 

Next we study how our results are changed when 
different sets of proton structure functions are cm- 
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Fig. 2. Rapidity distributions of the various channels for 
resolved photoproduction and of direct photoproduction in 
LO at y/s — 298 GeV and pt — 5 GeV evaluated with the 
LO sets by GRV and MT; see caption of Fig. 1 

ployed. Toward this end, we calculate d 2 a / dy dpT 
in the Born approximation at y/s = 298 GeV and 
y lab = —2 as a function of pt using in turn the proton 
structure function sets MT-B1 (PIS), MT-B1 (MS) 
[3] and MRSn-DO (MS), MRSn-D- (MS) [5], which 
are all of NLO. For the photon structure functions we 
choose the LO set by GRV [6] and for the QCD cou- 
pling a s we use the usual one-loop formula adopting the 
value of A from the respective proton structure func- 
tions. In Fig. 3, we show the cross sections normalized 
with respect to the calculation with the MT-SL set. 
We remark that only the MT-SL set should be com- 
bined with the LO hard-scattering cross sections be- 
cause it is a genuine LO fit; all the others correspond 
to NLO analyses of present data. 

We observe that the variation caused by replacing 
the LO set by the various NLO sets does not exceed 
30% for pt ft 2 GeV. A similar result was found in I 
for pp processes. We observe that the ratio of the two 
cross sections is always larger than one. From what 
was said above, it follows that the comparison of the 
NLO sets with the LO set is not particularly mean- 
ingful. However, it is instructive to compare the re- 
sults obtained with different NLO sets of proton struc- 
ture functions parametrizations. The corresponding 
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Fig. 3. Influence of typical NLO proton structure functions 
on the LO calculation of d 2 a/dydpT for resoved photopro- 
duction at y/s = 298 GeV and yi ab = —2. The curves 
are normalized to the LO calculation with the MT-SL set. 
In both numerator and denominator the GRV LO photon 
structure functions are used 
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Fig. 4. Influence of typical photon structure functions on 
the LO calculation of d 2 a/dydpT for resoved photoproduc- 
tion at y/s = 298 GeV and j/i a b = —2. The curves are 
normalized to the LO calculation with the GRV LO set. In 
both numerator and denominator the MT-SL proton struc- 
ture functions are used 



cross sections differ at most by 20%. The DIS and 
MS versions of MT-B1 differ by only 10% relatively 
to each other. We consider this as an upper bound on 
the scheme dependence, since the partial compensation 
that is expected to take place at NLO, due to scheme- 
dependent terms, is not in effect here. The MRSn 
DO and MRSn-D- sets fit also recent NMC [19] and 
CCFR [20] data, which extend to smaller values of x 
than previously tested. The DO and D— sets, available 
in the MS version, give cross sections with quite differ- 
ent behaviours at low px- At large px they give similar 
results and differ from the cross section obtained mak- 
ing use of the MT-B1 (MS) set by (10-15)%. When 
comparing with the pattern observed for pp reactions 
(see Fig. 2 in I), we must keep in mind that there the 
structure functions of both proton and antiproton are 
changed, whereas here the photon structure function is 
left unchanged. 

To study the influence of different photon structure 
functions, we calculate cPa/dydpT again in the Born 
approximation at y/s — 298 GeV and yi a b = — 2 as a 
function of px using in turn the NLO sets GRV (DIS7), 
GRV (MS) [6], ACFGP (massive) [8] and GS [7] and 
the LO parametrizations GRV [6], DG [18] and GS 
[7]. In all cases the MT-SL proton structure func- 
tion parametrization is used. The results are shown 
in Fig. 4. Except for low p T , the DG set gives cross 
sections which are 10% smaller than those for GRV 
(LO). The two steps in the DG curve are due to the 
fact that the DG parametrization is discontinuous at 
the charm and bottom thresholds, at px ~ 5 GeV and 



Pt ~ 14 GeV, respectively. The GRV structure func- 
tions are smooth at these thresholds. The GRV (DIS7) 
set gives results very close to the GRV (LO) results. 
This set involves a particular subtraction in the point- 
like part of the photon structure function, which is 
compensated by a change of the direct photoproduc- 
tion cross section. Since the latter is beyond the scope 
of this work, we shall not employ the GRV (DIS7) set 
in the following. We observe that the DIS7 and MS 
versions of the GRV parametrization differ apprecia- 
bly only at larger px, where direct photoproduction 
comes into play; at small pt (pt ~ 25 GeV) their dif- 
ference is below 10%. The cross section for the ACFGP 
choice of structure functions, which is also based on the 
MS scheme, is very different, especially in the lower-px 
range, where it is up to 50% smaller than the GRV 
results. This can be traced to the much weaker gluon 
component in the ACFGP set. At large px, where the 
influence of the gluon component diminishes, the ratio 
approaches the value 1 again. The results obtained for 
the LO and NLO GS photon structure function are not 
included in Fig. 4. The LO soft-gluon option yields val- 
ues for the ratio between 1 and 1.2 for px ~ 16 GeV. In 
the upper range of the px values considered in Fig. 4 
the ratio comes down to 0.9. The NLO set, which is 
based on the MS scheme, gives results in the ballpark 
of those obtained with the DG and ACFGP sets; for 
Pt ^ 20 GeV the GS (MS) values are somewhat smaller 
than the DG values, for px ~ 30 GeV they are in line 
with the ACFGP values. 
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To investigate the influence of the higher-order 
terms of the hard-scattering cross sections, we restrict 
ourselves to pure photoproduction, replacing the right- 
hand side of (2.1) by 6(1 — x). We calculate Sa(jp — ► 
hX)/dydpT at yfs = 256 GeV (corresponding to typ- 
ical 20 GeV photons colliding with 820 GeV protons) 
and j/ i a b = —2 as a function of px with and without the 
higher order terms and plot the ratio in Fig. 5. This 
ratio should change only little when the photon spec- 
trum (2.1) is taken into account. A realistic analysis 
including the WWA and also fragmentation will be re- 
ported in the next sub-section. At this point, we wish 
to separate the effect of the higher-order terms. To 
that end, we use the same structure functions, namely 
the MS versions of GRV and MT-B1, and the two- loop 
expression for a s in both numerator and denominator. 
The renormalization scale of a s is adjusted according 
to the proton structure functions, i.e., A = 194 MeV, 
which is very close to the value A = 200 GeV (Np = 4) 
implemented in the GRV set. We perform the compu- 
tation for the three choices p = M = Mf = px/2, px, 
and 2px- From Fig. 5 we see that the ratio is practically 
independent of px for px <L 10 GeV. For smaller pt it 
increases with decreasing pt at a rate that is largest 
for the smallest scale. Its value ranges between 1.5 and 
2.5, which demonstrates the importance of the NLO 
corrections to the parton-level cross sections. These 
results are very similar to the ones reported in I for 
pp reactions at y/s = 630 GeV. Note that, in con- 
trast to Figs. 3 and 4, the same structure functions 
are used in both numerator and denominator, so that 
variations due to changes of scheme or parametrization 
largely cancel out. We observe a strong scale depen- 
dence, which is mainly due to the denominator, as will 
become apparent below. 

At this point, a few comments on spurious higher- 
order effects are in order. In Fig. 5, we used NLO struc- 
ture functions also in connection with the NLO cor- 
rections to the hard-scattering cross sections. Strictly 
speaking, such a procedure introduces certain contribu- 
tions beyond next-to-leading order in the case of GRV. 
For the choice /i = M = Mf — px we have explicitly 
checked the influence of the modifications necessary to 
eliminate these spurious higher-order terms. At small 
Pt (pt ~ 5 GeV), the ratio is increased by some 10%, 
but already at px ~ 7 GeV the increase becomes in- 
significant, until it vanishes at px ~ 12 GeV. Then the 
ratio is reduced by an amount which reaches at most 
5% at the upper end of the px-range shown in Fig. 5. 
Similar results are obtained for single-hadron produc- 
tion: the only effect of the realistic fragmentation func- 
tions is to reduce the values of px at which the features 
previously described take place. Since at present there 
is still appreciable freedom concerning the details of 
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Fig. 5. NLO calculation in the MS scheme of d 2 a/dydpT 
for 7p — » j + X at \fs — 256 GeV and ?/i a b = —2 normal- 
ized to the corresponding calculation with the NLO terms in 
the hard-scattering cross section omitted. In both numera- 
tor and denominator the MS sets of the GRV and MT-B1 
structure functions are used. The scales fj,, M, and Mf are 
identified and set equal to pr/2, pt, and 2jjt 

the photon structure functions, we do not consider it 
necessary to worry about this complication, leaving it 
for a future, more refined analysis. 

The scale dependence of numerator and denomina- 
tor of the ratio plotted in Fig. 5 is explicitly investi- 
gated in Figs. 6a, b for px fixed at 10 and 20 GeV, re- 
spectively. The three scales /x, M, and Mf are set equal 
to £px and £ is varied between 0.2 and 5.0. Again, we 
choose the MS sets by GRV and MT, y/s = 256 GeV, 
and j/ i a b = —2. The LO cross section (dashed line) 
depends more strongly on the scale £ than the NLO 
result (solid line) . The more rapid decrease of the LO 
cross section for increasing £ explains the increase of 
the ratio plotted in Fig. 5. The NLO cross section is 
particularly insensitive to £ for £ > 0.5. This holds true 
independently of px, although the absolute cross sec- 
tion strongly varies with px- The scale variation of the 
NLO cross section is well bounded in contrast to the LO 
cross section, which decreases monotonically with in- 
creasing £. At this point, one should keep in mind that 
the Mf dependence in the NLO hard-scattering cross 
section is not compensated as long as 5-function-type 
fragmentation is assumed. If we do not vary Mf, which 
would be more appropriate in this case, we obtain the 
dot-dashed curves in Figs. 6a, b. This shows less £ de- 
pendence than the LO curves but more than the NLO 
curves with variable Mf. It seems that the factor- 
ization scale dependence in the NLO hard-scattering 
cross section is mostly cancelled between initial and fi- 
nal state. The pattern observed in Figs. 6a,b is very 
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Fig. 6a. Scale dependence of d 2 a/dydpT for 7p — > j + X 
at -v/s = 256 GeV, j/i a b = -2, and p T = 10 GeV. The solid 
and dot-dashed lines represent the NLO results in the MS 
scheme with variable and fixed Mf, respectively. The MS 
sets of the GRV and MT-B1 structure functions are used. 
The dashed lines correspond to the calculation with the 
NLO terms of the hard-scattering cross sections omitted 



similar to our findings regarding pp collisions in I. 



3.2 Single-hadron 
sections 



inclusive cross 



In this sub-section we shall present our results for the 
inclusive production of single charged hadrons and neu- 
tral pions. For the calculation of these cross sections 
we used the same fragmentation functions as in Sect. 2. 
We choose the MT-B1 and GRV structure functions, 
both in the MS scheme, except in Fig. 7, where we 
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Fig. 7. Ratio of direct over resolved contribution to the 
cross section for production of charged hadrons 

take the MT-SL and GRV (LO) structure functions. 
As in the previous sub-section, the parameter A in the 
coupling constant a s is set equal to the value required 
by the proton structure function parametrizations, i.e. 
194 MeV for NLO evaluations and 144 MeV for the LO 
ones. 

First we consider the ratio of direct over resolved 
production of charged hadrons h = (h + + h~ )/2 in or- 
der to see where they are equal in magnitude. We 
restrict ourselves to the LO calculation and we set 
V« = 298 GeV, yi ab = -2, and /j, = M = M f = p T . As 
in the previous sections, M indicates both, the factor- 
ization scales for the photon and the proton structure 
functions. The result is shown in Fig. 7. The cross- 
ing point, above which the direct production of single 
hadrons exceeds the resolved one, is at pr ~ 13 GeV. 
This value is considerably lower than the crossing point 
obtained in the case of LO jet-production cross section. 
We do not expect the position of these points to change 
drastically when NLO corrections are included. Fig. 7 
shows clearly that for pt > 5 GeV the direct produc- 
tion of single charged hadrons is not negligible anymore 
and our results can only apply to data with a visible 
photon remnant jet. 

In Fig. 8 we show d 2 a/dydp^ for inclusive single- 
charged hadron production as a function of px for 
V« = 298 GeV, y iah = -2, and \i = M = M f = p T . 
The upper curve is the NLO result and the lower curve 
the LO cross section with the same structure functions, 
MT-B1 and GRV in MS, and a s (/i 2 )- K factors rang- 
ing between 2.5 and 1.6 are obtained in the pt range 
considered here. They are smaller at hight pt than 
the factors observed in the case of J-function fragmen- 
tation (see Fig. 5), for the same choice of scales, but 



7 



10' 



10" 



ep -> h + anything 



> 10" 

CD 

O 



\ 

b 



10" 



10 r 



io 1 - 



: I I | I 


1 1 ' 1 = 
fi = M = M f = p T - 




Ml— Bl (MS) 


\\ 


GRV (MS) 


r \\ 
r Vs = 298 GeV 




: Ytab = -2 
i 


1 , , 



ep -> 7T + anything 
i | i i | 



Pt (GeV) 

Fig. 8. Inclusive charged-hadron production cross section 
for ^/s = 298 GeV and yiab = —2. The solid line correspond 
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level 



bigger at high px- Moreover, the inclusion of realistic 
fragmentation functions for the produced hadrons in- 
duces a decrease of the cross section, for increasing px, 
faster than the one observed in Fig. 1. 

Qualitatively similar results, shown in Fig. 9, are 
obtained in the case of single- n° production. However, 
the values of the cross sections are here lower than the 
ones obtained in the case of charged-hadron produc- 
tion, for the same choice of scales. It should be re- 
minded at this point that the results presented in I for 
single- 7T° production in pp collisions could fit the ex- 
isting data only for relatively high scales. Similar con- 
clusions were reached in [9] when the ^-fragmentation 
functions [11] were used. 

In connection with the results in Figs. 8 and 9, it 
is of interest to know the influence on the cross section 
of the gluon originating from the photon. According 
to the results shown in Fig. 1, we expect this gluon 
contribution to be more important at small px- In 
Fig. 10 we plot the pt dependence of the ratio 
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ah 



d 2 cr/dydpT - d 2 o / dy dp T \ (no g/j) 
d 2 <r/dydpT 



(2.7) 



Pt (GeV) 

Fig. 9. Same as in Fig. 8 for single-7r ( production 

where the label (no g/j) indicates that the gluon com- 
ponent in the photon has been put to zero. R g / 7 varies 
from 60% at p T = 2 GeV to 30% at p T = 10 GeV. This 
result is essentially independent of whether NLO cor- 
rections are included or not. 

The importance of the quark component of the pho- 
ton for increasing pt should become evident when we 
compare single-charged hadron production in ep and 
pp collisions at the same \fs. Since, at large px, the 
point-like part of the photon-quark component is much 
stronger than the proton-quark distribution, we expect 
the pp cross section to fall off more rapidly than the ep 
cross section, for increasing px . This is indeed the case, 
as can be seen in Fig. 11 where we compare ep — > hX 
via resolved photons with pp — > hX at y/s — 298 GeV 
and a CM rapidity at which both cross sections are 
maximal, ycM = —0.3. The two curves in Fig. 11 
are NLO (upper curves) and LO (lower curves) predic- 
tions, respectively. The values of the ep cross section 
have been multiplied by the factor 3000 which makes 
ep and pp results coincide up to px ~ 4 GeV. Already 
at px = 10 GeV, though, the two sets of curves deviate 
visibly from each other and at the upper end of the px- 
range shown in Fig. 11 the pp cross section falls short 
of the ep one by at least a factor two. Of course, in 
order to verify this effect experimentally we need pp 
and ep data at the same y/s and the isolation of the 
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Fig. 10. Fraction R g /~f of the charged-hadron production 
cross section at y/s — 298 GeV and yi a b = —2 due to the 
gluon originating from the photon. The solid line indicates 
the NLO result and the dashed line the LO one with struc- 
ture functions and a s at the NLO level 

resolved production in the ep process. With direct ep 
production included, the effect should be even stronger. 

Another issue to be considered before proceeding 
further is whether the inclusion of a fifth flavour intro- 
duces any appreciable deviation from the results pre- 
sented so far. We have verified that, in spite of the 
non-negligible contribution of the bottom-quark to the 
photon structure function, the increase in the cross sec- 
tions is at most 2%, for the highest px values here con- 
sidered, when we assume Nf — 5. 

Finally, in Fig. 12 we compare our predictions from 
Fig. 8 with recent data obtained by the HI Collabo- 
ration at HERA [1]. These data extend only up to 
p\ — 12 GeV 2 and have quite large experimental un- 
certainties for the largest values of p\. Since the abso- 
lute normalization of these data points is not known, 
we have adjusted them with an overall factor so as to 
have satisfactory agreement with the theoretical NLO 
results at p\ <~ 1.4-1.7 GeV 2 . We see that the p\ shape 
of the theoretical results is consistent with the data. 
Of course, for firmer conclusions, more accurate data 
which extend to higher values of p\ are needed. 

To obtain more information on the scale depen- 
dence of the inclusive cross section, we have calcu- 
lated it as a function of the scale factor £ defined in 
the previous sub-section for px — 3, 5 and 10 GeV and 
\/~s = 298 GeV, yi a b = —2. The results are shown in 
Fig. 13, where we compare the NLO predictions eval- 
uated with MT-B1 and GRV in MS and two-loop a s 
(solid lines) with the LO predictions evaluated with 




io 2 r i 
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Fig. 11. Comparison of the cross sections of charged- 
hadron production in pp collisions and via resolved photons 
in ep collisions at HERA energy and the same CM energy 
and rapidity, y/s = 298 GeV and ycM = —0.3. The upper 
curves indicate the full NLO results, while the lower ones 
give the LO results obtained with structure function and 
a B at the NLO level 

MT-SL and GRV (LO) and one-loop a s (dotted lines). 
To show the effect of the O (tvi?) terms in the hard- 
scattering cross section, we also plot the results ob- 
tained for the LO predictions when structure function 
and a s are kept at the NLO level. 

We see quite clearly that the NLO cross section has 
less scale dependence than the two LO ones. Except 
than for the case px — 3 GeV, the NLO cross section 
develops a plateau for £ = 0.6. The turn-over for small 
£ at pt — 10 GeV is presumably an artifact due to the 
use of LO fragmentation functions. Both LO curves 
show a monotonic dependence on £, which for £ > 1 
is steeper than the dependence of the NLO cross sec- 
tion. We also observe that the LO cross sections (dot- 
ted curves) equal the NLO ones for £ between 0.6 and 
0.9, whereas at larger values of £ the two sets of curves 
are quite well separated. These latter feature has also 
been observed in I for single-charged hadron produc- 
tion in pp collisions. 

A comparison of Fig. 13 for px = 10 GeV with 
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Fig. 13. Scale dependence of d 2 a/dydpT for the resolved 
contribution to ep — > foX at yfs = 298 GeV, yiab = —2, 
and three different values of pr, Pt = 3, 5, 10 GeV. The 
solid lines indicate the NLO result, the dotted lines the LO 
result, and the dashed lines the LO result with structure 
functions and a s at the NLO level 

scattering cross section do not change their shape in 
any significant way. 



Fig. 6a, where the scale dependence has been investi- 
gated in the case of <5-function fragmentation and real 
photoproduction at y/s = 256 GeV, yields the follow- 
ing observations. Apart from the absolute size, the 
LO results exhibit roughly the same qualitative fea- 
tures, which is attributed to a rather weak scale depen- 
dence of the LO fragmentation functions. The situa- 
tion changes when NLO corrections are included. The 
drop-off at small £ is stronger in Fig. 13 than in Fig. 6a. 
Finally, higher scales are needed to obtain equal values 
of LO and NLO results when realistic fragmentation 
functions are included. 

So far we have concentrated on the discussion of 
Pt distributions for fixed yi a b = —2, where cross sec- 
tions are maximal. In Fig. 14 we show the differential 
cross section d 2 a/dy dpr as a function of the rapidity 
y for yfs = 298 GeV, (j, = M = M f = p T , and two 
different values of p^: pt — 5 and 10 GeV. We observe 
that these distributions are quite symmetric around 
2/iab ~ 2 and that the NLO corrections to the hard- 



4 Summary and Conclusions 

We have presented next-to-leading-order predictions 
for inclusive single-charged hadron and single- n° pro- 
duction via resolved photoproduction at HERA. We 
have confronted these predictions with first experimen- 
tal data from the HI Collaboration and found that 
both are in good agreement as for the slope in pp- To 
allow meaningful tests of the QCD-improved parton 
model and to probe quantitatively the parton struc- 
ture of the photon, more accurate measurements are 
needed. 

We find that the LO and NLO calculations coincide 
when the renormalization and factorization scales are 
chosen to lie just below pr- We consider this agreement 
as accidental, since there is a substantial discrepancy 
for other choices of scales. As expected, the inclusion 
of NLO corrections leads to a reduction in the scale 
dependence of the cross sections. Therefore, within 



10 



ep -> h + anything 



CD 

a 



OH 



Oh 



10 



io J 



10 



10 



Vs = 298GeV 
fj. = M = M f = p T 



: 5GeV 




p T = lOGeV 



10" 

-6 -4-2 2 

yiab 

Fig. 14. Rapidity distribution of d 2 a/dydpT for the re- 
solved contribution to ep — > hX at ^/s = 298 GeV, pt = 
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at the NLO level 

the scale variation in NLO, our results can be consid- 
ered as an absolute prediction for resolved photopro- 
duction of single hadrons in ep collisions. This is par- 
ticularly important in the case of tt° production, since 
such events provide a background for prompt-photon 
emission. However, precisely in the case of neutral 
hadrons the fragmentation functions used in our cal- 
culation may not be completely reliable. 

A further limitation of our analysis is related to the 
use of LO fragmentation functions. Already now we ob- 
serve that the dependence on the fragmentation scale 
is reduced when the NLO terms of the hard-scattering 
cross section are included. We expect the results to be- 
come even more stable when also NLO fragmentation 
functions are utilized. 

Finally, in order to disentangle experimentally re- 
solved and direct photoproduction, it is necessary to 
select small-px hadrons and/or to require the detec- 
tion of the photon remnant jet. At present, this does 
not represent a problem, since the data are still being 
accumulated at relatively low p T . However, as soon 
as data at larger p T will become available, the direct 
process will have to be evaluated in NLO and matched 
to the calculation presented here. 
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